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Evaluating the relic density of dark matter is an interesting possibility to constrain the parameter 
space of new physics models. However, this calculation is affected by several sources of un- 
certainty. On the particle physics side, considerable progress has been made in the recent years 
concerning the calculation of the annihilation cross-section of dark matter, which is needed in this 
context. In particular, within the Minimal Supersymmetric Standard Model, the theoretical un- 
certainty has been reduced through the calculation of loop corrections. The present contribution 
gives an overview over the achievements that have been made in QCD corrections to neutralino 
pair annihilation. The numerical impact is illustrated for a few examples. 
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1. Motivation 

Among the most compelling evidence for physics beyond the Standard Model of particle 
physics is the presence of cold dark matter (CDM) in our universe. Assuming the standard cosmo- 
logical ACDM scenario, recent cosmological measurements constrain its relic density to 

0.1053 < Hcdm/j 2 < 0.1193 (1.1) 

at 95% confidence level. The given limit is based on the 7-year data from the WMAP mission in 
combination with other observations [PJ. 

For a given new physics scenario, the relic abundance of a dark matter candidate in our universe 
can be evaluated by solving the Boltzmann equation 

— = -3Hn- (a ann v) (n 2 - n 2 eq ) , (1.2) 

which describes the evolution of the number density n of the dark matter candidate. The first term 
on the right hand side conveys the dilution due to the expansion of the universe, H denoting the 
(time-dependent) Hubble parameter. The second term is related to annihilation and co-annihilations 
of dark matter. Here, n eq is the number density in thermal equilibrium. The total (co-)annihilation 
cross-section a ann , multiplied by the relative velocity v, has to be convolved with the velocity 
distribution to obtain the thermal averaged cross-section (a ann v). After integrating the Boltzmann 
equation, Eq. (1.2), the relic density of dark matter is given by £2cdm^ 2 = m x no/p c , m x and no 
being the mass and present number density of the dark matter particle, respectively, and p c the 
critical density of the universe. This cosmological constraint gives additional information with 
respect to mass limits from collider searches and precision measurements such as the rare decay 
b — > sy or the anomalous magnetic moment of the muon (g — 2)^. 

In the Minimal Supersymmetric Standard Model (MSSM), one of the most convincing dark 
matter candidates is the lightest neutralino, which is stable if it is the lightest superpartner and 
if /?-parity is conserved. Several public tools exist, that perform a numerical calculation of the 
neutralino relic density for a given set of parameters. In these packages, such as Dark SUSY |Q] or 
micrOMEGAs [^J, all processes that contribute to the annihilation cross-section are implemented, 
but mostly only as a pure leading order calculation. Radiative corrections are included only for 
very sensitive quantities, e.g., the bottom Yukawa coupling. However, it is well known that next- 
to-leading order corrections can have sizeable impacts on cross-sections, and may in consequence 
alter the prediction of the dark matter relic density. 

Recently, several results for corrections to the dark matter annihilation cross-section have been 
published. Here, we focus on corrections of &(cl s ) [Q, ||, g], which are most important in case of 
neutralino pair annihilation into quarks. For other processes, such as annihilation into gauge bosons 
or leptons, electroweak corrections can be important in the same manner. Details and a numerical 
analysis for this case can be found in Refs. [f7|]. Moreover, one-loop corrections due to the exchange 
of a gauge or Higgs boson in the initial state have been investigated in Ref. Finally, let us note 
that the prediction of the relic density is affected by further uncertainties. 

In the next Section, a brief overview over the annihilation of neutralinos into quarks and the 
associated corrections of order a s is presented. Numerical examples follow in Sec. ||. Finally, 
conclusions are given in Sec. 01. 
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Figure 1: Tree-level processes contributing to neutralino pair annihilation into quark-antiquark pairs. An- 
nihilation is possible through the exchange of a Higgs-boson (H? = h ,H°,A ), a Z-boson, or a squark 
(qi = qi,q 2 ). 



2. Neutralino annihilation into quark-antiquark pairs 

In the following, we focus on neutralino pair annihilation into quark-antiquark pairs. The 
relevant Feynman diagrams are depicted in Fig. [jj In the constrained MSSM (cMSSM) with unifi- 
cation conditions at the GUT scale, the annihilation into third generation quarks dominates in wide 
ranges of the WMAP-favoured parameter space [§J. At low values of the gaugino mass parameter 
mi/2, the annihilation proceeds efficiently through a resonance of the light scalar Higgs-boson h° 
into bb final states. In the so-called focus point region at high scalar masses mo, important final 
states are tt that occur through the resonance of the heavy scalar Higgs H . Finally, for large tan/3, 
the so-called "A-Funnel" region with resonant annihilation into bb through a pseudoscalar Higgs 
A opens. 

While in the cMSSM, mostly Higgs-resonances are relevant, the exchange of a Z-boson or a 
squark can be enhanced if (at least some of) the unification conditions are relaxed. In consequence, 
the higgsino fraction of the neutralino can be altered, leading to enhanced couplings to the Z-boson. 
This is the case, e.g., for scenarios with non-universal Higgs masses (NUHM) or non-universal 
gaugino masses. The exchange of squarks is preferred for light squark masses, which occurs in 
particular for a large trilinear coupling Aq. 

The annihilation cross-section related to the diagrams in Fig. [I] receives corrections of order 
a s due to one-loop diagrams with gluons and gluinos. These include self-energies for the quarks 
and squarks, corrections to the vertices involving quarks or squarks, and several box diagrams. At 
the same order in a s ., the emission of real gluons from one of the (s)quarks that are present in the 
tree-level diagrams has to be taken into account. For the full set of virtual and real contributions 
the reader is referred to Refs. 

All relevant amplitudes have been calculated analytically in Refs. [Q, |[ gj. After renormal- 
ization, the virtual part is ultraviolet-finite, the remaining infrared singularities cancel when taking 
into account the real gluon emission. The obtained results have been implemented in a numeri- 
cal package, which can serve as extension for the public relic density tools DarkSUSY |Q] and 
micrOMEGAs |gp. Note that we also take into account important effects due to the resumma- 
tion of the Yukawa coupling, in particular running quark masses and the resummation of terms 
proportional to tan/3 or the trilinear couplings A;,. Details on the calculation as well as concern- 
ing the renormalization procedure and the numerical implementation using the dipole subtraction 
formalism, can be found in Refs. [Q, ||, ^] and the references therein. 
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tan|3 m 1/2 (GeV) 

Figure 2: Impact of the next-to-leading order corrections on the annihilation cross-section for a typical 
scenario in the "A-Funnel" region (left) and the prediction of the neutralino relic density in the mQ-m l ^ 2 
plane for tan j3 = 54, Aq — 0, and > (right). The shown contours correspond to 0.094 < £2cdm^ 2 < 0.136. 

3. Numerical examples 

In the following, the numerical effect of the corrections discussed above is illustrated for a 
few examples. First, we consider the case of large tan/3 within the constrained MSSM. In the 
left panel of Fig. || we show the annihilation cross-section given by the diagram XiXi — ^ A — >• 
bb including different corrections for a typical reference point in the WMAP-preferred region of 
parameter space. The numerical values are normalized to the cross-section at leading order. The 
QCD corrections decrease the cross-section by more than 60%, which is mainly an effect of the 
running bottom-quark mass. Finite terms up to &(af) contribute by about 5%. 

The numerical impact of the SUSY-QCD corrections strongly increases with tan/3. While the 
effect is almost negligible for small tan /3 ~ 2, the cross-section is decreased by another 10% for 
tan /3 > 50, which is the case for the A-Funnel region. Again, the main impact comes from the 
resummation of the bottom Yukawa coupling, while finite terms are numerically less important [Q]. 

On the right panel of Fig. ||, we show the WMAP-favoured region of the mo-m l j 2 plane for 
fixed tan /3 = 54, Ao = 0, and /J. > and for different levels of correction. The relic density has been 
evaluated using DarkSUSY 4 . 1 which includes by default the QCD running quark mass and 
finite terms up to €?(aj) in the bottom Yukawa coupling, but no SUSY-QCD corrections. The 
latter have been implemented into the package in order to evaluate their numerical impact. For 
reference, we also show the prediction of the relic density including a pure tree-level calculation of 
the annihilation cross-section. As can be seen, the impact of the additional corrections is larger than 
the uncertainty of the WMAP measurement, especially for low mass parameters. Interestingly, the 
effect is reversed around the Higgs-pole at m\h ~ 800 GeV. This is explained by the fact, that the 
decay width of the Higgs-boson receives the same corrections as the annihilation cross-section [Q]. 

Let us now turn to the case where the annihilation proceeds mainly through the exchange of a 
Z-boson or a squark into a top-antitop final state. For a typical scenario within the NUHM scenario, 
we show in Fig. || the contributions of the different diagrams at the tree-level as well as the impact 
of the NLO corrections discussed in Sec. ^ on the cross-section. The shaded regions indicate the 
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Figure 3: The different contributions to the annihilation cross-section at tree-level (left) and effect of the 
discussed corrections (right) for a typical scenario with non-universal Higgs masses (left). The shaded area 
indicates the shape of the velocity distribution of the dark matter particle. 



velocity distribution needed for the thermal average. As can be seen, the Z- and squark-exchanges 
are linked by a strong destructive interference, so that the corrections to both subprocesses play 
an important role. At the same time, the contribution from Higgs exchange is negligible at the 
relevant energies. The total QCD and SUSY-QCD corrections increase the annihilation cross- 
section at tree-level by about 50% [0]. Note that the cross-section included in the public code 
micrOMEGAs 2 . 1 [^] takes into account parts of the NLO correction by using effective quark 
masses. 

Finally, in Fig. |] we show two examples of the preferred regions of parameter space, projected 
in the planes of the physical mass parameters m,A-\l and m t j-m«o, respectively. The correspond- 
ing high-scale scenarios feature non-universal Higgs masses and non-universal gaugino masses, 
respectively. The graphs have been obtained by scanning in m^-m^ and M\-mo planes, respec- 
tively, and projecting the contours into the planes of physical mass parameters. Here, the neu- 
tralino relic density has been evaluated using micrOMEGAs 2 . 1 ||3|], which has been linked to 
the numerical package discussed in Sec. || containing the full &(a s ) corrections to the annihilation 
cross-section. The three bands correspond to the levels of correction discussed above and shown in 
Fig-! 

As also for the A-Funnel region, the impact of the corrections is numerically more important 
than the uncertainty of the WMAP data. In our examples, we observe shifts of up to 50 GeV in 
the mass m A of the pseudoscalar Higgs boson, up to about 20 GeV for the mass m 9 o of the lightest 
neutralino, and up to about 200 GeV for the mass m~ tx of the lighter stop. Since the PLANCK 
satellite will deliver even more precise measurements in a very near future, the next-to-leading 
order corrections become even more significant. 



4. Conclusion 



In summary, we have presented a complete calculation of neutralino pair annihilation into 
third-generation quarks at next-to-leading order in a s . The result has been implemented in a nu- 
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Figure 4: Preferred regions in the niA-fi and mj -m^a planes for parameter variations around two typical 
scenarios within non-universal Higgs and non-universal gaugino masses, respectively. 



merical package, which allows to evaluate their impact on the cross-section. Via a link to public 
tools like DarkSUSY and micrOMEGAs the effect of the corrections on the prediction of the relic 
density has been studied numerically. We have shown that the discussed QCD corrections can have 
a sizeable numerical impact, in particular for scenarios that are inspired by minimal supergravity 
such as the constrained MSSM or models with non-universal Higgs or gaugino masses [Q, £0. 
The corrections should therefore be taken into account in the extraction of supersymmetric mass 
parameters from cosmological data or the analysis of the parameter space of supersymmetric theo- 
ries. Let us emphasize again, that both DarkSUSY and micrOMEGAs take into account a part of 
the presented corrections in terms of effective Yukawa couplings. 

The obtained results can easily be generalized to the case of the annihilation of two different 
neutralinos or charginos into quark pairs. Moreover, annihilation involving into light quarks can 
be included While such processes are negligible within "mSUGRA-like" scenarios, they can 
play an important role in more general setups such as the 19-parameter phenomenological MSSM 
(pMSSM), where the soft-breaking parameters are given at the electroweak scale. 

As already stated in the Introduction, electroweak corrections can also have a sizeable impact 
on the prediction of the neutralino relic density, in particular for annihilation into final states with 
gauge bosons and leptons [f7|]. Further uncertainties affecting the prediction of the relic density 
include uncertainties on the mass spectrum of the superpartners due to renormalization group run- 
ning [10] and loop-corrections to sparticle masses [11], assumptions on the cosmological model 
itself [12], and uncertainties on the effective degrees of freedom due to the QCD equation of state 
[0. 

Finally, let us note that the effect of the next-to-leading order corrections discussed in this 
contribution are independent from the other sources of uncertainty mentioned above. Moreover, 
higher-order corrections to the annihilation cross-section will become even more important when 
the PLANCK satellite will provide more precise cosmological data in a very near future. 

I would like to thank M. Klasen and K. Kovaffk for the fruitfull collaboration. This work is (in 
part) supported by the Landes-Exzellenzinitiative Hamburg. 
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